Reducing vortex density in superconductors using the ratchet effect 
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A serious obstacle that impedes the application 
of low and high temperature superconductor (SC) 
devices is the presence of trapped flux ||l|J^ . Flux 
lines or vortices are induced by fields as small 
as the Earth's magnetic field. Once present, vor- 
tices dissipate energy and generate internal noise, 
limiting the operation of numerous superconduct- 
ing devices [0,0. Methods used to overcome this 
difficulty include the pinning of vortices by the 
incorporation of impurities and defects |Q], the 
construction of flux dams ||5|, slots and holes ||^ 
and magnetic shields ||2|J^ which block the pene- 
tration of new flux lines in the bulk of the SC or 
reduce the magnetic fleld in the immediate vicin- 
ity of the superconducting device. Naturally, the 
most desirable would be to remove the vortices 
from the bulk of the SC. There is no known phe- 
nomenon, however, that could form the basis for 
such a process. Here we show that the applica- 
tion of an ac current to a SC that is patterned 
with an asymmetric pinning potential can induce 
vortex motion whose direction is determined only 
by the asymmetry of the pattern. The mecha- 
nism responsible for this phenomenon is the so 
called ratchet effect |@^5 and its working prin- 
ciple applies to both low and high temperature 
SCs. As a flrst step here we demonstrate that 
with an appropriate choice of the pinning poten- 
tial the ratchet eflect can be used to remove vor- 
tices from low temperature SCs in the parameter 
range required for various applications. 

Consider a type II superconductor film of the geometry 
shown in Fig. [|, placed in an external magnetic field H. 
The superconductor is patterned with a pinning poten- 
tial U{x, y) = U(x) which is periodic with period i along 
the X direction, has an asymmetric shape within one pe- 
riod, and is translationally invariant along the y direction 
of the sample. The simplest example of an asymmetric 
periodic potential, obtained for example by varying the 
sample thickness, is the asymmetric sawtooth potential, 
shown in Fig. Rb. In the presence of a current with den- 
sity J flowing along the y axis the vortices move with the 
velocity 



(fL + fv 



fu)/r/, 



(1) 



pointing in the direction of the external magnetic field 
H, fu — — ^x is the force generated by the periodic 
potential, fw is the repulsive vortex-vortex interaction, 
$0 = 2.07 X 10~^ G cm^ is the flux quantum, r] is the vis- 
cous drag coefficient, and d is the length of the vortices 
(i.e. the thickness of the sample). 

When a dc current flows along the positive y direction, 
the Lorentz force moves the vortices along the positive 
X direction with velocity u+. Reversing the current re- 
verses the direction of the vortex velocity, but its mag- 
nitude, |w-|, due to the asymmetry of the potential, is 
different from Vj^. For the sawtooth potential shown in 
Fig. |l|b the vortex velocity is higher when the vortex is 
driven to the right, than when it is driven to the left 
(f+ > l^-l). As a consequence the application of an ac 
current (which is the consecutive application of direct 
and reverse currents with density J) results in a net ve- 
locity V = {v+ + v-)/2 to the right in Fig. nb. This net 
velocity induced by the combination of an asymmetric 
potential and an ac driving force is called ratchet veloc- 
ity l|3"0 • The ratchet velocity for the case of low vortex 
density (when vortex-vortex interactions are neglected) 
can be calculated analytically. Denoting the period of 
the ac current by T, the ratchet velocity of the vortices 
in the T ^ oo limit is given by the expression 
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if /l < /i 

if /i < /l < /2 

if /2 < /l, 
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where fL = (J x h)<I>od/c is the Lorentz force moving the 
vortices transverse to the current, h is the unit vector 



where /i — tJJ jix and /2 = AC// £2 are the magnitudes 
of the forces generated by the ratchet potential on the 
facets of length (.\ and £2, respectively (see Fig. ^3), tJJ 
is the energy difference between the maximum and the 
minimum of the potential, and /l = |fL| = J'^od/c. 

Since for high magnetic fields vortex-vortex interac- 
tions play an important role, we have performed molecu- 
lar dynamics simulations to determine the ratchet veloc- 
ity for a collection of vortices. As Fig. || demonstrates, 
we find that for low vortex densities the numerical results 
follow closely the analytical prediction (^ , and the mag- 
nitude of the ratchet velocity decreases with increasing 
vortex density. The vortex densities used in the simula- 
tions correspond to an internal magnetic field of about 
0.7, 35, and 70 G, covering a wide range of magnetic 
fields. A key question for applications is if the ratchet 



velocity (g) is large enough to induce observable vortex 
motion at experimentally relevant time scales. To ad- 
dress this issue in Fig. g we plotted v for Nb, a typical 
low temperature SC used in a wide range of devices, for 
which the potential U{x) is induced by thickness varia- 
tions of the SC. The details of the model are described in 
the caption of Fig. g. As the figure indicates, the maxi- 
mum ratchet velocity (5.2 m/s) is high enough to move 
a vortex across the typical few micrometer wide sample 
Q in a few microseconds. Furthermore, increasing the 
vortex density by two orders of magnitude decreases the 
vortex velocity only by a factor of three. 

Next we discuss a potentially rather useful application 
of the ratchet effect by demonstrating that it can be used 
to drive vortices out of a SC. Consider a SC film that is 
patterned with two arrays of the ratchet potential ori- 
ented in opposite directions, as shown in Fig. ^. During 
the application of the ac current the asymmetry of the 
potential in the right half moves the vortices in that re- 
gion to the right, while vortices in the left half move 
to the left. Thus the vortices drift towards the closest 
edge of the sample, decreasing the vortex density in the 
bulk of the SC. We performed numerical simulations to 
quantitatively characterize this effect, the details and the 
parameters being described in the caption of Fig. 0. In 
Fig. |3|b we summarize the effectiveness of vortex removal 
by plotting the reduced vortex density inside the SC as a 
function of the Lorentz force /l and the period T of the 
current. One can see that there is a well defined region 
where the vortex density drops to zero inside the SC, in- 
dicating that the vortices are completely removed from 
the bulk of the SC. Outside this region we observe either 
a partial removal of the vortices or the ac current has no 
effect on the vortex density. 

The (/L,r) diagram shown in Fig. ^3 has three ma- 
jor regimes separated by two boundaries. The Ti — 
'^V f _rf +f.^ i-yj+i )] phase boundary (here we assume 
d/2 < £2 and /in (a;) is defined in Fig. y) provides the 
time needed to move the vortex all the way up on the 
ti long facet of the ratchet potential at the edge of 
the SC, i.e. to remove the vortex from the SC. When 
T < Ti the vortices cannot exit the SC. The T2 = 

2'? ( /.-[/'+/oa.o] + f,-[fJTfi-l+d/2)] ) Pliase boundary 
(where /edge is also defined in Fig. ^) is the time needed 
for a vortex to enter from the edge of the SC past the 
first potential maxima. Thus, when T < T2 the vor- 
tices cannot overcome the edge of the potential barrier. 
These phase boundaries, calculated for non-interacting 
vortices, effectively determine the vortex density in the 
three phases. Vortex removal is most effective in regime 
I, where the vortices cannot move past the first potential 
barrier when they try to enter the SC, but they get past 
the barriers opposing their exit from the SC. Thus the 
vortices are swept out of the SC by the ratchet effect, and 
no vortex can reenter, leading to a vortex density p — 0. 
Indeed, we find that the numerical simulations indicate 
complete vortex removal in the majority of this phase (see 



the contour lines in Fig. gb). An exception is the finger 
structure near the crossing of the Ti and T2 boundaries. 
For fields and periods within the first finger (lowest in 
Fig. Ip) the vortex follows a periodic orbit inside a sin- 
gle potential well flG]. The subsequent fingers represent 
stable period orbits between two, three, or more wells, 
respectively. Since the vortices cannot escape from these 
orbits, they remain trapped inside the SC, increasing the 
vortex density within the fingers in the phase diagram. 
Fig. ||b shows the analytically calculated envelopes of the 
regions where such trapping occurs. An important fea- 
ture of the finger structure is that stable periodic orbits, 
which prevent vortex removal, do not exist above the line 
Ttip = fh f fit _f \ connecting the finger tips. In regime 
II vortices can enter the SC, but the ratchet effect is still 
sweeping them out, thus here we expect partial removal 
of the vortices, the final vortex density inside the SC be- 
ing determined by the balance of vortex nucleation rate 
at the edge of the sample (which depends on the surface 
properties of the SC) and the ratchet velocity moving 
them out. In regime III the vortices cannot leave the SC 
and new vortices cannot enter the system, thus the ini- 
tial density inside the SC is unchanged throughout this 
phase (p = po)- 

Since the forces fin{x) and /edge depend on H, the po- 
sition of the phase boundaries Ti and T2 also depends on 
the external magnetic field. In particular, there exists a 
critical field H* , such that for H > H* regime I, where 
vortex removal is complete, disappears, but regime II 
with partial vortex removal does survive. We find that 
for Nb films of geometry described in Fig. we have 
H* « lOG. However, since H* is a consequence of the 
geometric barrier, its value can be modified by changing 
the aspect ratio of the film. Furthermore, for supercon- 
ductors with elliptic cross section the geometric barrier 
can be eliminated [ll|, thus phase I with complete vor- 
tex removal could be extended to high magnetic fields as 
well. 

Vortex removal is important for numerous SC applica- 
tions and can improve the functioning of several devices. 
An immediate application of the proposed method would 
be improving the operation of superconducting quantum 
interference devices (SQUIDs), used as sensors in a wide 
assortment of scientific instruments ||,|l^,|l^. A long- 
standing issue in the performance of SQUIDs is 1// noise 
[g,|2[ , arising from the activated hopping of trapped vor- 
tices ||l|. Reducing the vortex density in these super- 
conductors is expected to extend the operation regime of 
these devices to lower frequencies. 

Although over the past few years several applications 
of the ratchet effect have been proposed, such as sepa- 
rating particles | |l^ , |l5|l, designing molecular motors [^ , 
smoothing surfaces [|17f|, or rectifying the phase across a 
SQUID p8| , p^ , our proposal solves a well known acute 
problem of condensed matter physics, by removing vor- 
tices from a SC. In contrast with most previous applica- 
tions, which require the presence of thermal noise, our 



model is completely deterministic. Indeed, in Nb the 
variations in the pinning potential is A[/ w 25 eV, which 
is more than 10'' times larger than kBT « 0.8 meV at 
Tc = 9.26 K, thus rendering thermal fluctuations irrel- 
evant. On the practical side, a particularly attractive 
feature of the proposed method is that it docs not re- 
quire sophisticated material processing to make it work: 
First, it requires standard few-micron scale patterning 
techniques (the micrometer tooth size was chosen so that 
a few teeth fit on a typical SQUID, but larger feature size 
will also function if the period T is increased proportion- 
ally). Second, the application of an ac current with ap- 
propriate period and intensity is rather easy to achieve. 
For applications where an ac current is not desired, the 
vortices can be flushed out before the normal operation 
of the device. On the other hand, if the superconduct- 
ing device is driven by an ac current (e.g. rf SQUIDs, 
ac magnets, or wires carrying ac current), the elimina- 
tion of the vortices will take place continuously during 
the operation of the device. The analytically predicted 
phase boundaries, whose position is determined by the 
geometry of the patterning, provide a useful tool for de- 
signing the appropriate patterning to obtain the lowest 
possible vortex density for current and frequency ranges 
desired for specific applications. Finally, although here 
we limited ourselves to low temperature SCs, the working 
principle of the ratchet effect applies to high temperature 
superconductors as well. 
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FIG. 1. Patterning with an asymmetric potential of a SC. 
(a) Schematic illustration of a SC in the presence of an ex- 
ternal magnetic field _H". A dc current with density J fiowing 
along the y direction (indicated by the large arrow) induces 
a Lorentz force /l that moves the vortex in the x direction. 
The SC is patterned with a pinning potential U{x, y) = U(x), 
whose shape is shown in the lower panel. The potential is 
periodic and asymmetric along the x direction, and is trans- 
lationally invariant along y. (b) The pinning potential U{x) 
along the SC cross section. The solid arrows indicate the vor- 
tex velocity «+ («_) induced by a direct -I- J (reversed —J) 
current. The average, v = {v+ + V-)/2, is the ratchet veloc- 
ity of the vortex, obtained when an ac current is applied, (c) 
The parameters characterizing a single tooth of the asymmet- 
ric potential. 
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FIG. 2. Ratchet velocity of the vortices as a function of the 
amplitude of the driving force /l • The thick solid line corre- 
sponds to the analytical result (bl) for a single vortex line. The 
symbols are the result of the numerical simulations for multi- 
ple vortices. The simulations were done using the model de- 
veloped by Nori and collaborators [pC)|-p2| , assuming that the 
rigid vortices are pointlike object moving in the x-y plane. At 
time zero the vortices are positioned randomly in the SC with 
a density po and they move with velocity given by (fi]). The 
vortex-vortex interaction between two vortices at position r^ 

and Vj is modeled using [gj fw ~ g^°^i Ki{{ri — Vj)/\)vij, 
where iij = (r^ — r j ) / 1 r^ — Tj | . Here the modified Bessel func- 
tion Ki is cut off beyond the distance r = 25A, where A is 
the penetration depth (for Nb A = 45nm at T = 0). The 
force / generated by the sawtooth pinning potential shown 
in Fig. hi is equal to /i when kt < x < k£ + £i, and /2 
when k£ + ii < X < {k + 1)1, where fc = 0, 1, ..., iV - 1. We 
choose £i = 20A = O.dfim, ^2 = 5A = 0.225Atm, £ = I1+I2 and 
A^ = 10, giving for the total width of the sample w — 11.25^m. 
Its length (along the y direction) is set to 12^m. The sample 
has periodic boundary conditions in both the x and y direc- 
tion. The Lorentz force due to the ac current is equal to +/l 
for T/2 time, and -/l for T/2 using T = O.S/is. We con- 
sidered the simplest case, in which the potential is induced 
by thickness variations of a Nb superconductor thin film of 
thickness d, i.e. the SC thickness, d + h{x), changes along 
the X direction, following a sawtooth pattern. The pinning 
energy acting on the vortices is given by U{x) = {d + h{x))eo, 
where eo is the line energy of the vortex per unit length. 
Thus the magnitudes of the forces acting on the vortices are 
/i = eoAh/li and /2 — eoA/1/^2 for the two facets of the Ah 
high teeth (shown in Fig. hic), and we choose A^ = £2- For 
Nb we have eo = 1-7 x 10~ ^N, the viscosity per unit length 
is »?o = 7 X 10"®Ns/m^, yielding -q = rjod = 1.4 x 10"^^Ns/m 
for a d = 2OOOA thick film. The total number of vortices in 
the simulation were n = 5(o), n = 250(l>), and n — 500(o) 
corresponding to ~ 0.7G(o), 35G(l>), and 70G(o) magnetic 
field in the sample. 



(a) 



V 



U(x) 



kkKkmaa^aa 



V 



X 



■w 



w 




p/ f\ 

1.0 



I 



0.8 
0.6 
0.4 
0.2 
0.0 



0.0 1.0 2.0 3.0 4.0 5.0 6.0 

/l(x10^'N) 



FIG. 3. Removing vortices from a SC using an asymmet- 
ric potential, (a) The potential necessary to remove vortices 
from the SC. Using the simulation method described in Fig. bl 
we investigated a system consisting of A'' = 5 teeth oriented to 
the left and the same number oriented to the right, as shown 
in the figure, the parameters of each tooth being identical to 
that described in Fig. hp. To mimic the pressure generated by 
the external magnetic field, which acts to push vortices into 
the SC, on the two sides we attached two reservoirs, that have 
a constant vortex density po at all times. Thus, vortices can 
leave the SC for the reservoir, or new vortices can enter from 
the reservoir. In thin SC films, due to the Meissner current, 
there is a geometrical barrier that acts to trap the vortices 
inside the SC pj. Since most applications of SC's involve 
thin films, we included in the simulations this geometrical 
barrier, that creates a force fin{x) — — ^x/\/w^ — x^ for 

-w + d/2<x <w- d/2, and /cdgo = 2eo - ^y/'iw/d-l 
for X > w — d/2, and —/edge for x < ~w + d/2. Thus the geo- 
metrical barrier opposes the entry of the vortices at the edge 
of the SC, but once they move inside, it moves them towards 
the center of the SC. For successful vortex removal the ratchet 
effect has to be strong enough to move the vortices against 
fin{x). (b) The {fh,T) diagram describing the effectiveness 
of the ratchet effect as a function of the parameters charac- 
terizing the driving current, /l. The color code corresponds 
to the the relative vortex density p/po, where po is the initial 
vortex density corresponding to H = IG and p is the final 
vortex density after the application of the ac current. As the 
color code indicates, there is a region where vortex removal is 
complete, the vortex density being equal to zero. The dashed 
lines correspond to the Ti and T2 boundaries, that are cal- 
culated analytically (see text) and separate the three main 
regimes: I: complete vortex removal in the majority of the 
regime, p = 0; II: partial vortex removal, < p < po; and 
III: no change in the vortex density, p = pa. The thin white 
lines denote the boundaries of the regions where vortex trap- 
ping due to periodic orbits occurs. These boundaries correctly 
reflect the structure of the fingers, but slightly deviate from 
the results of the numerical simulation, because the analytical 
calculation assumed an array of identical teeth. 



